Hyperfine structure in hydrogen and helium ion 
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QED theory of the hyperfine splitting of the Is and 2s state in hydrogen isotopes and helium-3 
ion is considered. We develop an accurate theory of a specific difference 8-Bhfs(2s) — -E"hfs(1s)- We 
take into account fourth order corrections and nuclear structure effects. The theoretical prediction 
is now of a higher accuracy than the experiment is. The study of the difference provides the most 
accurate test (on a level of a part in 10 s ) of the QED theory of Is HFS up to date. The theory 
agrees with most of the experimental data. 
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The hyperfine structure (HFS) intervals of the ground 
state in a number of neutral atoms and singly charged 
ions can be measured with a high accuracy. However, 
theory even in the case of the simplest of them (such as 
hydrogen isotopes and the helium-3 ion) is essentially af- 
fected by nuclear structure effects which contribute from 
30 to 200 ppm and cannot be calculated accurately. In 
contrast the Is HFS interval in muonium is calculated 
with an uncertainty of about 0.1 ppm and can be used 
to accurately test the bound state Quantum Electrody- 
namics (QED). However, the muonium calculations in- 
volves precision values of the fundamental contants (a 
and m^/rrie) and it would be important to test the QED 
calculations for the HFS without interferring with such 
a problem. 

Study of a specific difference 



D 2 i =8£hfs(2s)-£ H fs(1s) 
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provides us with an opportunity to make a test of the 
QED theory on a level of accuracy essentially better than 
1 ppm §. Such a high accuracy is possible because of 
an essential cancelation of nuclear contributions. We re- 
port here on new results for the difference in Eq. (Q). 
We complete calculations of the fourth order corrections 
and find nuclear structure contributions which remain 
after cancelation of the leading effects. Some of the cor- 
rections obtained here are bigger than the experimental 
uncertainty and must be taken into account. 

Our results are found to be in a fair agreement with 
most experimental data on hydrogen, deuterium and 
helium-3. We present a significant improvement of the 
theory for D21 and demonstrate that the comparison of 
theory and experiment for helium tests presently 

the QED theory of Is and 2s HFS on the highest level, 
namely one part in 10 s . That superceeds the muonium 
HFS by an order of magnitude. 

The hyperfine splitting of an ns state in a hydrogen- 
like atom with a nuclear mass M and a nuclear spin / 
can be presented in the form 
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Here Ry is the Rydberg constant, c is the speed of light, 
h is the Planck constant, /is is Bohr's magneton and m 
is the electron mass. The nuclear magnetic moment /i in 
our notation can be negative (if its direction is opposite 
to the nuclear spin) and the Fermi energy Ep, related to 
an energy splitting between the atomic state with total 
angular moment F = 1+1/2 and 7—1/2, can be negative 
as well. The QED correction for the HFS interval in the 
ground state is (see Ref. j| for references) 
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where a e is the electrons anomalous magnetic moment. 
A comparison of the QED calculations with experimental 
values is summarized in Table 3. To compute theoretical 
values we use fundamental and auxiliary constants from 
Refs. |,|. The QED expression above does not take into 
account any recoil effects. Recoil contributions involve 
high momentum transfer fl6f| and are essentially affected 
by the nuclear structure. In Table | we also present data 
for the 2s state, the theoretical expression for which is 
similar to Eq. ([I]) but some coefficients are different (see 
below). 

One can see that the Is hyperfine structure has been 
measured very accurately but any test of the QED cal- 
culations is limited by an essential contribution related 
to the nuclear structure which cannot be calculated pre- 
cisely. In fact the uncertainty of the nuclear-structure 
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contribution is at least 20% in hydrogen p8[ , and for 
deuterium the accuracy is not better |1S|| , In the case of 
tritium and helium-3 ion no results on the contribution 
of the nuclear effects has been obtained to the best of our 
knowledge. Thus, the pure QED theory is incomplete be- 
cause of lack of the nuclear-structure contributions and 
a comparison of the QED theory with the experiment in 
Table j] demonstrates how much it is incomplete. Our fi- 
nal target is a comparison of Is and 2s HFS intervals and 
for this reason we do not try to correct the QED theory 
for the nuclear effects. Contrary, we compare the pure 
QED calculation and experimental data to "measure" the 
nuclear contribution. 

For comparison we presented in Table | a theoretical 
result on Is muonium HFS jl^] , which contains the re- 
coil contributions and even small non-QED terms. Muo- 
nium, being a pure leptonic atom, is free of the nuclear- 
structure problem, however, the accuracy of any theo- 
retical calculation is limited to 10~ 7 by the uncertainty 
of experimental values for parameters needed to calcu- 
late the Fermi energy in Eq. (||). Those are the muon 
magnetic moment and the fine structure constant. Be- 
low we demonstrate that combining data for the Is and 
2s hyperfine structure in hydrogen and 3 He + we can go 
far beyond 1 ppm level ||] and hence develop a precision 
test of the QED theory for HFS compatible with the one 
related to muonium HFS. 

The theory for the specific difference D21 in Eq. (|l|) 
up to the third order in units of the Fermi energy was 
developed some time ago po|-^2| 
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The nuclear-structure corrections essentially shift the 
HFS value from its QED prediction. Three major nuclear 
effects contribute to the difference 



A^(ls) = E e *» s (ls) - £§™(ls) 
in Table |[ Namely they are: 



(G) 



• the nuclear charge and magnetic moment distribu- 
tion (that is the biggest effect in the case of hydro- 
gen); 

• a nuclear polarizability contribution (that is the 
biggest effect in the case of deuterium); 

• nuclear recoil contributions of order (Za)(m/M)EF 
and higher. 



There is also a correction to the Lamb shift caused by 
the nuclear structure 



A<^(bs) = -AZafm^Rl , 



(7) 



where Re is the nuclear electric charge radius and rel- 
ativistic units in which h = c = 1 are used. When the 
contributions to HFS (|J) and the Lamb shift (fjj) are de- 
termined, one can try to obtain a correction for differ- 
ence £>2i- That is possible because most of the nuclear- 
structure corrections do not depend on the details of the 
atomic structure. Both, the leading contributions to the 
HFS and the Lamb shift are of a special factorized form 



A£(Nucl) = A(Nucl) x \^ n i(r = 0)1 



(8) 



The energy shift is the product of the nuclear-structure 
parameter A(Nucl) and the value of the wave function 



|^(r = 0)| 2 = ^^o 



(9) 



The leading correction to the difference in Eq. (Q) must 
therefore vanish. The non- vanishing contributions can be 
expressed in terms of some effective <5-like potentials 



T/(Nucl) = A(Nucl) • 8(r) 



(10) 



The coefficient A(Nucl) can be for various nuclear contri- 
butions calculated (see e.g. Eq. (Q)) or determined from 
a comparison of experiment and a pure QED theory (se e 
e.g. Eq. (|)). The result is of the form (cf. Refs. @,||!) 
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where 1 + £ = R 2 M /R 2 E is a ratio of quadratic magnetic 
and electric nuclear radii. We obtain the nuclear struc- 
ture contribution to the Is HFS interval from comparison 
in Eq. (^|) and conservatively estimate the uncertainty as 
10%. The Lamb shift contribution is taken from Eq. (Q). 

The fourth order corrections to D 2 \ have been inten- 
sively studied for the last three years. The logarithmic 
corrections in order a 2 (Za) 2 and a(Za) 2 (m/M) were 
calculated in Ref. @ (^ HtH), the (Za) 3 (m/M) terms 
are found below (cf. |24| ). The only term calculated pre- 
viously is the relativistic term of the order [Zaf'Ep p6| . 

Partial results on the a(Za) 3 contributions were found 
in Ref. Q. They are related to effective non- relativistic 
potentials which lead to logarithmic contributions for the 
Is state HFS (25). The terms in the same order should 
also appear from potentials which contain some deriva- 
tives. A complete result on the self energy contribution 
was calculated after a suggestion by us in Ref. p7[ . We 
report here the completion of the evaluation of the vac- 
uum polarization effects. We derive an exact result for 
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HFS of the 2s state and a contribution to D21 is found 
via a comparison with the previously obtained result for 
the Is state ||. 

The fourth order contibutions are finally found to be 
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where 



C S£ (Z = 1) = 2.07(25) , C SE (Z = 2) - 2.01(19) 
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The partial results for the constants Cse and Cyp that 
are obtained in Ref. pj contain some misprints. Being 
corrected, the partial results (Cse — 2.5 and Cyp — 
0.83) are found to be close to the complete results above. 
That confirms an intuitive assumption that the potentials 
with derivatives lead to relatively small contributions. 
Smallness of terms with derivatives is important for our 
estimation of uncertainties of the nuclear-structure cor- 
rections. 

Let us discuss the uncertainty of the QED expression. 
The first two terms in Eq. ( |l2] ) are found in the logarith- 
mic approximation and we estimate the next-to-leading 
terms by a half-value of the leading contribution. How- 
ever, in the case of the third term in Eq. ( 12 ) the situation 
is more complicated. First of all, the (Za)[in/M)Ep cor- 
rections to the Is HFS contain a nuclear-structure depen- 
dence presented by ln(mi?£). Since we have not included 
them into the QED expression (H), they are effectively 
taken into account as a part of AE^pg . That means that 
an essential part of the (Za) 3 (m/M)Ep contribution into 
D21 is effectively included into Z?2i(Nucl) via Eq. (11). 
However, there are some contributions with loop momen- 
tum of about one electron mass and below which does not 
depend on the nuclear structure. They can be enhanced 
because of a relatively big magnetic moment (compared 
to the Dirac value) and we estimate the uncertainty of 
the last term in Eq. (gj) as (jjl/ fi B )(Za) 3 E F (cf. Eq. (|)). 

All contributions to the difference D21 in hydrogen, 
deuterium and helium-3 ion are summarized in Table ||. 
Parameter £ is known very badly, but it is not expected 
to be much larger than unity and hopefully the £-term 
is essentially below the uncertainties related to theory 
and experiment and thus may be excluded from further 
considerations. 



An essential improvement of the theory is achieved. In 
previous papers related to third-order QED corrections 
pO| , pH the uncertainty was not spelled out. We found 
here a number of corrections exceeding the experimental 
uncertainty. We state that after the examination pre- 
sented here the theoretical predictions (Table |l|) 



£>2i(theor) = D^f (QED) + D 2 \> (QED) + L> 21 (Nucl) 

are more accurate than the experiment. Five accurate 
measurements performed on three atomic systems are 
compared with our calculation in Table III. Four ex- 



(4), 



perimental results are in fair agreement with our the- 
ory, but a recent result for hydrogen jll| shows a 1.8cr 
discrepancy. The most important comparison is related 
to 3 He + : the 2s HFS was measured most accurately 
and its value is also the most sensitive to higher order 
corrections (because of larger Z and larger nuclear con- 
tributions). Because of a fair agreement of our theory 
with the helium experiment we expect that in the case of 
hydrogen the discrepancy is related to a problem on the 
experimental side. 

We consider comparison of theory and experiment for 
the difference D21 as a test of a calculation of a state- 
dependent part of corrections to E®p S ( ns) and hence 



we present in the last column in Table III a standard 



deviation a with respect to the Fermi energy Ep, i. e. 
to a value directly related to the Is HFS. That compar- 
ison demonstrates that study of D21 in helium ion pro- 
vides a more accurate test of QED than the study of the 
muonium HFS (a/Ep ~ 0.1 ppm) and indeed of HFS 
in hydrogen and other atoms with a structured nucleus. 
The uncertainty for D21 in hydrogen and deuterium is 
determined experimentally, while in the case of helium a 
value of a contains an essential contribution from theory 
as well. 

Most of so-called QED tests involve in part some other 
problems such as 

• verification of nuclear models and calculations of 
nuclear effects and hadronic contributions; 

• tests of consistency of data for fundamental con- 
stants (such as muon magnetic moment) or effec- 
tive parameters (such as the proton charge radius) 
related to completely different experiments. 

The D21 theory is free of all these problems. No con- 
stants are really involved: an effective value of £^pg'(ls) 
related to the nuclear effects arises from HFS theory Its 
contribution being relatively small is under control as 
well as other nuclear contributions. 

It is important to mention that presently there 
are three crucial higher-order QED contributions to 
hydrogenic energy levels: radiative recoil of the 
order a(Za) 2 (m/M)Ep, pure recoil of the order 
(Za) 3 (m/M)Ep and two-loop effects of the order 
a 2 (Za) 6 m. The difference D21 is sensitive to all of them 
and a progress in its calculation will therefore contribute 
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into progress in theory of the hydrogen Lamb shift, muo- 
nium hyperfine structure and positronium energy levels. 
The most accurate measurement on this difference is re- 
lated to an 25-years old experiment on helium ions and 
we can hope that some experimental progress to improve 
the most precise test of QED theory for the hyperfine 
structure is possible. 
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TABLE I. Hypefine structure in light atoms. Hydrogen re- 
sult for Is j?J is an average value over the most accurate data 
||. The difference between experiment and QED theory is 
denoted by AE. It is related to the nuclear contribution. 



Atom, 
state 


pexp 
-^HPS 

[kHz] 


pQED 
HFS 

[kHz] 


n 3 AE/E F 
[ppm] 


H, Is 


1420 405.751768(1), jR 


1420 452 


- 33 


D, Is 


327384.352 522(2), IT 


327 339 


138 


T, Is 
3 He+, Is 


1516 701.470 773(8), |0 


1516 760 


- 38 


- 8 665 649.867(10), ]§ 


- 8 667 569 


222 


H, 2s 


177 556.785(29), 111] 


177 562.7 


- 33 


H, 2s 


177 556.860(50), H 




- 32 


D, 2s 


40 924.439(20), (if 


40 918.81 


137 


3 He+, 2s 


- 1083 354.981(9)7H 


- 1 083 594.7 


221 


3 He+, 2s 


- 1083 354.99(20), R 




221 


Mu, Is 


4 463 302.78(5), |J 


4 463 302.91(56) 


0.0(1) 



TABLE II. Various contributions to L>2i(theor). Theoret- 
ical predictions depend on parameter £ = R 2 M /R 2 E — 1. 



Value 


H 


D 


3 He+ 


D^i (QED) [kHz] 
D^i (QED) [kHz] 
£>2i(Nucl) [kHz] 


48.937 
0.018(3) 

- 0.002 

- 10" 4 c 


11.305 6 
0.004 3(5) 
0.002 6(2) 
- 10" 4 C 


- 1 189.262 

- 1.137(53) 
0.331(36) 
+ 0.009 C 


D 2 i(theor) [kHz] 


48.953(3) 

- io- 4 c 


11.312 5(5) 
- IO" 4 c 


- 1 190.067(63) 
+ 0.009 C 



TABLE III. Precision tests of QED theory for D 2 i. The 
final standard deviation a includes contributions from both: 
theory and experiment. AD21 stands for the difference of ex- 
periment and theory. References for the D21 are presented 
for experiment for the 2s HFS. We put here ( — 0. 



Atom 


D21 (exp) 


D21 (theor) 


AD21/0- a/E F 




[kHz] 


[kHz] 


[ppm] 



H 
H 
D 

3 He 4 
3 He 4 



48.53(23), [11 
49.13(40), [12 
11.16(16), @ 
1 189.979(71) ji 
-1190.1(16), ff 



48.953(3) 

11.312 5(5) 
1 190.067(63) 



- 1.8 
0.4 

- 1.0 
0.9 
0.0 



0.16 
0.28 
0.49 
0.01 
0.18 
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